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ABSTRACT 

We present a new method to compute the cosmic emissivity Ej, and background intensity Jj,. Our 
method is based entirely on data from quasar absorption-line studies, namely, the comoving density 
of HI and the mean metallicity and dust-to-gas ratio in damped Lya galaxies. These observations, 
when combined with models of cosmic chemical evolution, are sufficient to determine the comoving 
rate of star formation as a function of redshift. From this, we compute Ej, and Jj, using stellar 
population synthesis models. Our method includes a self-consistent treatment of the absorption and 
reradiation of starlight by dust. In all of our calculations, the near-UV emissivity declines rapidly 
between z Ki 1 and z = 0, in agreement with estimates from the Canada-France Redshift Survey. The 
background intensity is consistent with a wide variety of observational limits and with a tentative 
detection at far-IR wavelengths. 

Subject headings: cosmology: diffuse radiation — galaxies: evolution — quasars: absorption lines 



1. INTRODUCTION 

The mean emissivity of the universe 5^ and the mean 
intensity of background radiation are important cos- 
mological probes. The former is defined here as the 
power radiated per unit frequency per unit comoving 
volume, while the latter is defined as the power received 
per unit frequency per unit area of detector per unit solid 
angle of sky. These quantities, of course, are not inde- 
pendent; is given by an integral of 5^ over redshift. 
It is likely that 5^ and are dominated at near-UV, 
optical, and near-IR wavelengths by the direct radia- 
tion from stars and at far-IR wavelengths by reradiated 
starlight from the dust within galaxies. Thus, they con- 
tain potentially valuable information about the global 
history of star formation. Observationally, 5^ and 
have proven to be elusive; until recently, it was only 
possible to estimate 5^ at z ^ 1 and to place weak 
constraints on J^. Theoretically, they have also proven 
to be elusive. Most calculations of 5^ and are based 



on the properties of present-day galaxies and some as- 
sumed evolution in the past, often specified by several 
free parameters (see Lonsdale 1995 for a review). 

In this Letter, we introduce a new method to compute 
5^ and J^; in essence, we predict the "emission history" 
of the universe from its "absorption history." From the 
absorption lines in the spectra of distant quasars, it is 
possible in principle to determine the global rates of gas 
consumption, metal production, and hence star forma- 
tion in galaxies at redshifts up to z 4. The focus 
here is on the damped Lya galaxies, which contain most 
of the cool, neutral gas in the universe and appear to 
be the progenitors of present-day galaxies (Lanzetta, 
Wolfe, & Turnshek 1995). We have already combined 
absorption-line data with models of cosmic chemical evo- 
lution to compute the comoving densities of stars and 
dust in damped Lya galaxies (Pei & Fall 1995). Here, we 
employ stellar population synthesis models to compute 
5^ and J^, including the absorption and reradiation of 
starlight by dust. In this first illustration of the method. 
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we have deliberately kept the analysis as simple as possi- 
ble in order to highlight the main ideas and assumptions. 
Moreover, the existing data on damped Lya galaxies are 
too sparse to warrant a more elaborate analysis. 

2. MODELS 

Our first task is to relate the cosmic emissivity and 
background intensity, Ej, and J^, to the comoving den- 
sities of stars and dust, and (expressed here in 
units of the present critical density, pc = ZHq/SttG). 
We ignore other sources of radiation (primarily active 
galactic nuclei) and exclude photons that are absorbed 
in the galaxies in which they were produced. Thus, for 
the stellar part of the emissivity, we write 

= (1 - A,)£l, (1) 

where A,, is the mean fraction of absorbed photons, and 
£g^, is the emissivity before absorption. The latter is 
given by 

CW = Pc^*dtX(t-OY^^' (2) 

where S„(At) is the power radiated per unit frequency 
per unit initial mass by a generation of stars with an 
age At, R is the returned fraction, and the dot denotes 
differentiation with respect to cosmic time. We com- 
pute S„ and R using the latest version of the Bruzual- 
Charlot models (described by Chariot, Worthey, & Bres- 
san 1996). In most of the calculations reported here, 
the stellar initial mass function (IMF) is assumed to be 
a power law, ^(m) cx m~(^"'"^), with x = 1.5 and upper 
and lower cutoffs at lOOM© and O.IMq. We further 
assume that all H ionizing photons are absorbed in the 
local interstellar medium, that 68% of them are con- 
verted to Lya photons (the fraction appropriate for case 
B recombination in gas at lO'* K), that all of these are 
absorbed by dust (as a consequence of resonant scatter- 
ing by HI), and that the remaining energy is radiated 
uniformly in wavelength between 3000 and 7000 A (a 
range that includes most of the relevant emission lines). 
This treatment of ionizing radiation is consistent with 
far-UV observations of starburst galaxies (Leitherer et 
al. 1995). 

We relate the mean fraction of photons absorbed by 
dust to the other properties of galaxies as follows. For 
simplicity, we assume that the stars and dust have the 
same spatial distributions within galaxies (so that the 
source function is constant along any ray), and we ignore 
the influence of scattering on absorption (which should 
be a good approximation for average quantities such as 
j4^). In this case, the fraction of photons absorbed along 



a ray that intercepts a galaxy with an optical depth 
is given by 

a(T,) = l-T-i[l-exp(-T,)]. (3) 

We now define ^(T^)dT^ to be the fraction of such rays 
with optical depths between and + dr^ when all 
galaxies and all positions and directions within them are 
considered. This enables us to express the mean fraction 
of photons absorbed by dust in the form 

Jo dT„T„(p(T„) 

The optical depth and HI column density N along a ray 
are related by = K,,kmmHN , where is the opacity 
(i.e., mass absorption coefficient), km is the dust-to-HI 
mass ratio, and mjj is the mass of an H atom. Thus, 
for a single value of km, the (true) distributions of 
and TV" are related by (j)(T,,) cx f(N). For consistency 
with the models of chemical evolution, we adopt the 
gamma distribution, f(N) = (f^./N)ex.p(-N/N^.) (Pei 
& Fall 1995). This implies 4i(t^) oc exp(— t^/th,^) 
with Th,^ = Ki^kmrnnN^,. Equations (3) and (4) and the 
relations 0,; = fcrnOfji and Ohi = (S'KGmH /•icHo)f*N^, 
then give 

A^ = l- T-J ln(l + Th.^), (5) 

Th.^ = (3cffo/87rG/H.)K^O<j. (6) 

If the stars and dust have different spatial distributions 
or there is a dispersion in the dust-to-gas ratio, these 
expressions should be regarded as approximations. We 
have checked that our results are not sensitive to the 
particular form of A,, . 

The starlight absorbed by dust is reradiated ther- 
mally. We approximate the corresponding emissivity by 
that of a single blackbody: 

£d. = 47rpeO<jK.S,(T<,). (7) 

The effective temperature of the dust is then deter- 
mined by the condition of energy balance 

47rpeO<j / diyK,[B,{Td) - B,(Tcmb)] = / diyA,£°„ 
Jo Jo (8) 

where Tcmb = 2.73(1 + z)K is the temperature of the 
cosmic microwave background radiation. Equation (7) 
should be a good approximation at wavelengths beyond 
the peak [A^ax ~ 150(T(j/20 K)~^ /^ni], where the emis- 
sion is likely to be dominated by large amounts of cool 
dust, but a poor approximation at shorter wavelengths, 
where the emission is likely to be dominated by small 
amounts of warm dust (e.g., small, transiently heated 
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grains). We return to this point later. For the opti- 
cal properties of the dust, we adopt the Draine & Lee 
(1984) model but with the proportions of graphite and 
silicates adjusted so as to reproduce the mean UV ex- 
tinction curve in the Large Magellanic Cloud (Pei 1992). 
In this model, the opacity at long wavelengths has the 
form Ki, oc v'^, and the left-hand side of equation (8) is 
proportional to 0(j(T^ — Tcf^g). Given 0, and 0,;, we 
can now compute the emissivity Ej, = Egi, + Edv from 
the equations above, and the corresponding background 
intensity at z = from 

dt 



J. 



dz£i 



dz 



(9) 



Here, we have neglected absorption between the sources 
of radiation and the observer. This is appropriate be- 
cause ionizing photons are assumed to be absorbed lo- 
cally and because relatively few non-ionizing photons are 
absorbed by the dust in intervening galaxies. 

The comoving densities of stars and gas in galaxies, 
and Og, and the mean metallicity in the interstellar 
medium Z, including dust, are governed by the equa- 
tions of cosmic chemical evolution: 



9.gZ - yn, = (Zf - Z){lf. 



(10) 
(11) 



Here, y is the IMF-averaged yield, and the source terms 
on the right represent the inflow or outflow of gas with 
metallicity Zf at a rate ^f. For purposes of illustration, 
we assume that just over half of the heavy elements are 
locked up in dust grains and that the ionized and molec- 
ular components of the interstellar medium are negligi- 
ble, i.e., 0(j = 0.55.^Og and 0^ = 1.30hi- The con- 
stancy of the dust-to-metals ratio can also be expressed 
in the form k{z)lZ{z) = fc(0)/^(0) = 0.8/^®, where k is 
the ratio of Tg, the extinction optical depth in the rest- 
frame B band, to TV/IO^^ cm'^ (see Table 2 of Pei 1992). 
This agrees to within a factor of two with the dust and 
metal content of present-day galaxies and damped Lya 
galaxies at z = 2.2 (Pei, Fall, & Bechtold 1991; Pettini 
et al. 1994). Given Ohi and some assumptions about 
Oji and Zf, it is now possible to solve equations (10) and 
(11) for and 0,;. We adopt the solutions presented 
by Pei & Fall (1995). These are of three types: a closed- 
box model (Oji = 0), a model with inflow of metal-free 
gas (Oji = +jyOs, Zf = 0), and a model with outflow 
of metal-enriched gas (Oji = — j^O^, Zf = Z). The ad- 
justable parameters are the initial comoving density of 
gas in galaxies Ogoo and the relative inflow or outflow 
rate v. The models include self-consistent corrections 
for the damped Lya galaxies that are missing from op- 



tically selected samples as a result of the obscuration of 
background quasars (Fall & Pei 1993]). The models were 
designed to reproduce the observed comoving density of 
HI in damped Lya galaxies at < z ^ 4 (Lanzetta et 
al. 1995, Storrie-Lombardi, McMahon, & Irwin 1996). 
They also reproduce the observed mean metallicity Z 
0.1.^0 at z = 2.2 (Pettini et al. 1994) and are consistent 
with the average properties of present-day galaxies. The 
results presented here, unless otherwise noted, are based 
on models with Ogoo = 4 x 10~^h~^, v = 0.5, and h = 
0.5, go = 0.5, A = (with h = ffo/100 km s"! Mpc'^). 
For comparison, we also display results from the analo- 
gous models without dust. 

3. RESULTS 

Figure 1 shows the evolution of 5^ at rest-frame wave- 
lengths of 2800 A and 1.0 /im. The former is dominated 
by the light from young, massive stars and is thus nearly 
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Fig. 1. — Cosmic emissivity 5^ as a function of red- 
shift z at rest-frame wavelengths of 2800 A (top) and 
1.0 /im (hoiiom). The solid curves are from the closed- 
box (C), inflow (I), and outflow (O) models with dust, 
while dashed curves are from the analogous models with- 
out dust. The data points with error bars are estimates 
from local samples and the Canada- France Redshift Sur- 
vey (Lilly et al. 1996). 



a direct indicator of current star formation, while the 
latter includes the light from a wide range of stellar 
types and thus reflects a combination of past and cur- 
rent star formation. The solid curves in Figure 1 rep- 
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resent the closed-box (C), inflow (I), and outflow (O) 
models with dust, while the dashed curves represent the 
analogous models without dust. The data points with 
error bars are estimates by Lilly et al. (1996) from lo- 
cal samples of galaxies (z = 0) and from the Canada- 
France Redshift Survey (z = 0.3, 0.6, and 0.9). These 
include corrections for incompleteness at faint magni- 
tudes. The observed emissivity at 2800 A decreases by 
an order of magnitude between z = 0.9 and z = 0, in- 
dicating a similar decrease in the comoving rate of star 
formation. Evidently, the slopes of the predicted and 
observed £i,(z) relations agree to within the uncertain- 
ties, both at 2800 A and 1.0 /im. The amplitudes match 
for the closed-box and outflow models with dust but not 
for the inflow model with dust. However, the predicted 
amplitudes are somewhat sensitive to the adopted IMF 
in the population synthesis models. If the IMF slope 
were increased to x = 2.0 or the lower cutoff were de- 
creased to 0.01 M0, the predicted amplitudes would be 
reduced by factors of 2 — 3. With this freedom to adjust 
the IMF, the emissivities in the inflow model with dust 
can be brought into rough agreement with the observed 
emissivities. We count this as a success because models 
with very different histories of star formation cannot be 
made consistent with the observations for any choice of 
the IMF. For example, models with O^oo ^ 1 x 10~^/i~^ 
or Ogoo ^ 8 X 10~^h~^ have emissivities that decline too 
slowly or too rapidly. 

Figure 2 shows vj,, as a function of wavelength. 
Again, the solid curves represent the closed-box (C), in- 
flow (I), and outflow (O) models with dust, while the 
dashed curves represent the analogous models without 
dust. The symbols with arrows and the hatched line 
represent the observational limits described in the cap- 
tion. Evidently, our models are consistent with these 
constraints, over four and a half decades in wavelength. 
The near-UV background is higher in the models with 
dust because they have more star formation at low red- 
shifts than the models without dust (see Figure 1 and 
Pei & Fall 1995). The open circles with vertical bars 
in Figure 2 represent the tentative detection of an ex- 
tragalactic far-IR background by Puget et al. (1996). 
This result, derived from COBE/FIRAS data, is uncer- 
tain because it depends critically on the removal of fore- 
ground emission by interplanetary and interstellar dust. 
Our models with dust are consistent with this detection. 
The effective temperature of the dust remains near 
20 K at z ^ 1 and then decreases to T,; 15 K at z = 0. 
This produces a peak in vj,, at A 240 /im. The valley 
at A 60 /im is an artifact of our assumption that the 
spectrum of the dust emission is that of single blackbody 
[see equation (7)]. We have experimented with more re- 
alistic, two-temperature models and find that the valley 
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Fig. 2. — Background intensity at z = times 
frequency as a function of wavelength A. The solid 
curves are from the closed-box (C), inflow (I), and out- 
flow (O) models with dust, while the dashed curves are 
from the analogous models without dust. The squares 
with arrows are lower limits derived from galaxy counts 
at 2000 A (Milliard et al. 1992), 3600-9000 A (Tyson 
1995), and 2.2 /im (Cowie et al. 1994). The triangles 
with arrows are upper limits or possible detections at 
1600 A (Paresce 1990) and 4000 A (Mattila 1990). The 
unadorned arrows are the residuals of COBE/DIRBE 
data after preliminary removal of foreground emission at 
1.25-240 /im (Hauser 1996). The hatched line is a con- 
servative upper limit based on COBE/FIRAS data at 
500-5000 fim (twice the limit from Mather et al. 1994; 
see Hauser 1996). The open circles with vertical bars are 
a tentative detection derived from COBE/FIRAS data 
at 170-1260 /im (Puget et al. 1996). 

can easily be filled in by emission from small amounts 
of warm dust with little effect on the emission at longer 
wavelengths. This indicates that the far-IR background 
is a robust feature of our models. The amplitude of 
depends on the IMF in the same way as the amplitude of 
£„; in particular, the background would be weaker at all 
wavelengths if the proportion of massive stars were re- 
duced. Moreover, the values of Ogoo that are consistent 
with the data on are essentially the same as those 
that are consistent with the data on 5^. 

4. DISCUSSION 

We have computed the cosmic emissivity 5^ and back- 
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ground intensity Jj, with input only from absorption-line 
studies of damped Lya galaxies. These objects have 
N ^ 10^*^ cm~^ (by definition) and are the probable 
sites of most star formation in the universe. They appear 
to be the progenitors of present-day galaxies, although 
many of their properties remain to be determined. In 
particular, we do not yet know the sizes and morpholo- 
gies of the damped Lya galaxies. It is possible that most 
of those at low redshifts are disks, while most of those at 
high redshifts are spheroids. We emphasize that the re- 
sults presented here are not affected by such issues. The 
reason for this is that all of the quantities required in our 
analysis can be computed directly from statistics of the 
absorption along random lines of sight [e.g., Ohi is given 
by the integral of N f(N) over TV"]. The primary restric- 
tion on our results is that they do not include any galax- 
ies that consumed their HI before z 4, the highest 



redshift probed systematically by quasar absorption-line 
studies. Since our calculations include several approxi- 
mations and idealizations, we regard them as illustrative 
rather than definitive. Nevertheless, we find reasonable 
agreement with all of the available data on 5^ and 
without any fine-tuning of parameters. In particular, 
our models reproduce the rapid decline in the near-UV 
emissivity between z 1 and z = reported by Lilly et 
al. (1996). They are also consistent with a wide variety 
of observational limits on the extragalactic background 
and a tentative detection at far-IR wavelengths reported 
by Puget et al. (1996). It is therefore possible that this 
background is produced mainly or entirely by galaxies 
at z ^ 4. 

We thank M.G. Hauser, N. Panagia, and the referee 
for helpful comments. 



REFERENCES 



Chariot, S., Worthey, C, & Bressan, A. 1996, ApJ, 457, 
625 

Cowie, L.L., Gardner, J. P., Hu, E.M., Songaila, A., Ho- 
dapp, K.-W., & Wainscoat, R.J. 1994, ApJ, 434, 114 
Draine, B.T., & Lee, H.M. 1984, ApJ, 285, 89 
Fall, S.M., & Pei, Y.C. 1993, ApJ, 402, 479 
Hauser, M.G. 1996, in Unveiling the Cosmic Infrared 
Background, ed. E. Dwek (Woodbury: AIP Press), 
11 

Lanzetta, K.M., Wolfe, A.M., & Turnshek, D.A. 1995, 
ApJ, 440, 435 

Leitherer, C, Ferguson, H.C., Heckman, T.M., & Lowen- 
thal, J.D. 1995, ApJ, 454, L19 

Lilly, S.J., Le Fevre, O., Hammer, F., & Crampton, D. 
1996, ApJ, 460, LI 

Lonsdale, C.J. 1995, in Extragalactic Background Radi- 
ation, ed. D. Calzetti, M. Livio, & P. Madau (Cam- 
bridge University Press), 145 

Mather, J.C. et al. 1994, ApJ, 420, 439 

Mattila, K. 1990, in lAU Symposium 139, The Galac- 



tic and Extragalactic Background Radiation, ed. S. 
Bowyer & C. Leinert (Dordrecht: Kluwer), 257 
Milliard, B., Donas, J., Laget, M., Armand, C, & 

Vuillemin, A. 1992, A&A, 257, 24 
Paresce, F. 1990, in lAU Symposium 139, The Galac- 
tic and Extragalactic Background Radiation, ed. S. 
Bowyer & C. Leinert (Dordrecht: Kluwer), 307 
Pei, Y.C. 1992, ApJ, 395, 130 
Pei, Y.C, & Fall, S.M. 1995, ApJ, 454, 69 
Pei, Y.C, Fall, S.M., & Bechtold, J. 1991, ApJ, 378, 6 
Pettini, M., Smith, L.J., Hunstead, R.W., & King, D.L. 

1994, ApJ, 426, 79 
Puget, J.-L., Abergel, A., Bernard, J. -P., Boulanger, F., 
Burton, W.B., Desert, F.-X., & Hartmann, D. 1996, 
A&A, 308, L5 
Storrie-Lombardi, L.J., McMahon, R.G., & Irwin, M.J. 

1996, MNRAS, submitted 
Tyson, J. A. 1995, in Extragalactic Background Radia- 
tion, ed. D. Calzetti, M. Livio, & P. Madau (Cam- 
bridge University Press), 103 



5 



